Collective roughening of elastic lines with hard core interaction 
in a disordered environment 
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We investigate by exact optimization methods the roughening of two and three-dimensional sys- 
tems of elastic lines with point disorder and hard-core repulsion with open boundary conditions. In 
2d we find logarithmic behavior whereas in 3d simple random walk -like behavior. The line 'forests' 
become Eisymptotically completely entangled as the system height is increased at fixed line density 
due to increasing line wandering. 
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I. INTRODUCTION 

The paradigmatic model for the mixed or Shubnikov 
phase of a super-conductor in a magnetic field is a sys- 
tem of interacting elastic lines, representing the magnetic 
vortex lines threading the superconductor along the field 
direction P] ■ The number of lines per unit area perpen- 
dicular to the field axis, also called the line density p, 
is given by the magnetic field strength and also defines 
the typical line-to-line distance. The interaction between 
the lines is repulsive and long ranged but exponentially 
cut off beyond a screening length, A. As long as this 
screening length is larger than the typical line distance 
a = Xjp^l'^ (i.e. A > a) at low temperatures the lines 
arrange in a periodic pattern, the Abrikosov flux line lat- 
tice. In the presence of quenched disorder such as it is 
produced by crystal lattice defects or impurities this long 
range order is lost and for weak disorder or low magnetic 
field strengths transformed into quasi- long-range order, 
the Bragg glass phase For strong disorder the ar- 
rangement of lines is highly irregular and the system is 
in a state that is called the vortex glass Q . 

In this paper we want to study the case when the 
screening length A is much smaller than the typical line 
distance A ^ a). This corresponds to high magnetic 
field or strong disorder and small values of A. From a 
theoretical point this limiting case is particularly inter- 
esting since it corresponds to a situation in which the sys- 
tem cannot be described by an elastic theory any more, 
which is based on the assumption of weak disorder and 
low line density. Perturbative approaches that start from 
a ground state with long range order and assume small 
fluctuations turn out to be inappropriate in this strongly 
disordered situation and one has to rely on numerical 
calculations to study the low temperature properties. 

Here we study the extreme limit of a system of lines 
in a disordered environment with hard core repulsion. 
Although it is inspired by strongly disordered supercon- 
ductors in a magnetic field it describes the more general 
situation of a dense system of one-dimensional (line-like) 
objects that compete for the energetically most favorite 
locations in the matrix they live in but repel each other 
via a hard core exclusion principle - for instance long 



polymers stretching from one side of a very inhomoge- 
neous sample to the opposite side. We use a simplified 
(directed) polymer model and study the disorder aver- 
aged ground state properties of this system. To cal- 
culate the exact ground states we use an optimization 
algorithm, the minimum-cost-flow-algorithm, which de- 
termines the exact optimal solution in polynomial time 

El. 



We focus on the roughness properties of the lines (i.e 
the typical transverse fluctuations of the lines on their 
way from the top to the bottom of the system) and will 
show that in two dimensions the steric repulsion between 
the lines is sufficient to produce the same scaling behav- 
ior of the roughness as predicted by the elastic theories 
for long range interactions, weak disorder and low line- 
densities. In particular in 2d the hard core repulsion 
leads to collective rearrangements of the lines that yield 
a roughness that increases logarithmically with system 
size, also called super- rough behavior "qj, Q]- In three 
space dimensions the situation is totally different: Elas- 
tic theories predict that the roughness increases with 
the square-root of the logarithm of the system size 
whereas we find for the case we consider that the lines 
become more or less transparent for each other and can 
wander transversally from one side of the sample to the 
other. Steric repulsion alone is thus not sufficient to re- 
strict the transverse fluctuations of a line system in 3d. 



The paper is organized as follows: In the next section 
we will introduce the model and specify the methods by 
which we compute and analyze its ground states. In sec- 
tion 3 and 4 we will present our results for the roughness 
of the 2d and 3d system, respectively, and present the 
corresponding finite size scaling forms. In section 5 we 
consider various generalization of the disorder, in par- 
ticular weak and anisotropic randomness, and derive the 
scaling laws that govern the crossover behavior to the 
universal scaling forms reported in the previous sections. 
Section 6 contains a discussion and an outlook. 
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II. MODEL 

We consider the following model of a system of inter- 
acting lines in a two- or three-dimensional disordered en- 
vironment: The lines live on the bonds of a simple cubic 
lattice with a lateral width L and a longitudinal height 
H (i.e. M = LxHotM = LxLxH lattices sites in 
2d or 3d, respectively) with free boundary conditions in 
all directions. The lines start at the bottom plane and 
end at the top plane, if necessary the entrance and exit 
points can be fixed. The number N of lines threading 
the sample is fixed by a prescribed density p — N / L and 
p = N jl?' in two and three dimensions respectively. We 
restrict ourselves to hard-core interactions between the 
lines, which means that their configuration is specified 
by bond-variables rii G {0,1}, = 1 indicating that a 
line segment occupies a bond with index j, and ti^ = 
indicating that no line segment occupies this bond. 

We model the disordered environment by assigning 
a random (potential) energy ei G [0, 1] (uniformly dis- 
tributed) to each bond i, such that the total energy of 
the line configuration is given by 
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FIG. 1: Line identification schemes in 2 and 3 dimensions. 
In (a) and (b) the bold and dotted lines denote the choice 
identification of two lines, (c): an example of the case where 
lines are identified randomly. 



upon a bond configuration is not unique. Since we want 
to calculate the roughness of lines we need to determine 
the individual lines, for which we use a local rule. In 2d 
the line identification is unambiguous if we simply require 
that the lines cannot cross, see Fig. la. The same rule is 
applied in 3d when the line segments ending in a single 
site are within the same plane (see Fig. lb), otherwise 
the choice is random (see Fig. Ic). 



In section 5 the distribution of is modified from the 
presented above in order to extend our analysis also to 
the cases of weak and anisotropic disorder. 

In a continuum limit the system of interacting elastic 
lines is described by the following Hamiltonian: 
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dz 



+ Vr[n{z),z] 



+ J2Vinth{z)~r,{z)]}, (2) 

where ri{z) denotes the transversal coordinate at lon- 
gitudinal height z of the i-th flux line. Our numer- 
ical model corresponds to the case where the interac- 
tions T4nt[ri(z) — rj(z)] are hard core repulsive and the S- 
correlated disorder potential Vr[ri{z), z] has to be strong 
compared to the elastic energy that is proportional to 7. 

At low temperature the line configurations will be dom- 
inated by the disorder and thermal fiuctuations are negli- 
gible. Therefore we restrict ourselves to zero temperature 
and focus on the ground state of the Hamiltonian 
Computing the ground state now corresponds to finding 
N non-overlapping paths traversing the network from top 
to bottom. One has to minimize the total energy of the 
whole set of the paths and not of each path individually 
(already the two-line problem is actually non-separable 
Pll. This task can be achieved by applying Dijkstra's 
shortest path algorithm successively on a residual graph 

m 

Although the lines cannot occupy the same bond of the 
lattice they may touch in isolated points as exemplified 
in Fig. 1. This means that the line identification based 



III. ROUGHENING IN 2D 

The main quantity we are interested in is the disorder 
averaged line roughness. The mean square displacement 
of a single line with index i in one sample is defined as 



w. 



(3) 



with r"i = i^~^Ef=l^"(■^) for n = 1,2. The mean 
square displacement of all lines in one sample is = 

X^i^i ''^i ^^'^ f^"3 roughness w is defined as the square 
root of the disorder averaged sample mean square dis- 
placement w = \J {w'^). 

The roughness of a one line system {N — 1) scales as 
w ~ H'' in the limit of infinite transverse system size. In 
2d the value of the roughness exponent is (J = 2/3 jj], 
whereas in 3d it is close to ^ = 5/8 0, ^| • In the case of 
a non-vanishing line density one expects to observe this 
single line behavior as long as the transverse fiuctuations 
of the individual lines are smaller than the average line- 
to-line distance a, which is given by the line density p — 
1/a in 2d. This means that we expect w ~ H'' for H <C 

Once the transverse fluctuations of the individual lines 
have reached the size of average line-to-line distance one 
expects a collective behavior of the lines that restricts 
the individual line roughness due to the presence of the 
others. If the line system behaves like an elastic medium 
the roughness in the collective regime is expected to be- 
have like w ^ \nL Hence, for fixed line density p we 
expect the following scaling form 



w ^ a\n{L) ■ g2d{H/{a\nLy/^) 



(4) 
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FIG. 2: Scaling plot of the roughness in 2d according to the 
finite size scaling form I^J- The line density is p = 0.05, the 
roughness exponent is = 2/3. 



FIG. 4: Scaling plot for the roughness in 2d according to the 
scaling form JSJ. The system size is fixed ai L — 256. 



o p=0.1 y = 0.46 + 0.96 ln(x) 
□ p=0.2 y = 0.81 +0.431n(x) 
o p=0.3 y = 0.68 + 0.29 ln(x) 
A p=0.4 y = 0.63 + 0.20 ln(x) 
V p=0.5 y = 0.49 + 0. 17 ln(x) 






o □ — — ■ 

□ 




□ IX * ^ 

1 1 ^=^=^=^ 



10 100 

L 



FIG. 3: Saturation roughness as a function of the system 
width with fits of the form y = A + B\n{L). 



eral values of the line-density p. The data can be fitted 
to the form y = A + Bln{L) again reflecting the col- 
lective super-rough scaling of the roughness in 2d. The 
slope B of the data sets decreases as p increases while 
the constant part A does not vary much. This leads to 
the strengthening of finite size effects with increasing p. 

The crossover from single-line to multi-line scaling 
takes place when H'^ ~ a, where a is the average line- 
to-line distance a = 1/p in 2d. For fixed but large lateral 
system size L the scaling form Q predicts 

w^a- hdiH/a'/'^) (5) 

where the scaling function g2d{x) has the asymptotic be- 
havior g2d{x) ^ x'^ for a; <C 1 and g2d{x) ^ const, for 
X ^ 1. Fig. 21 shows the corresponding data collapse 
of the roughness data that we computed. Note that the 
height has been also rescaled with a factor of 1 — 1/a to 
account for the limit p — > 1. 



where g2d{x) is a scaling function with the asymptotic 
behavior g2d{x) oc x'' for a; <C 1, corresponding to the 
single-line behavior, and g2d{x) = const, for a; ^ 1, cor- 
responding to the collective regime. Note that we have 
assumed that the line density enters this form only via a 
rescaling of the lateral length scales. 

We computed the ground states of a large number of 
disorder realizations (up to 10'^) for various values of L, 
H and p and produced scaling plots for fixed line density 
p according to the suggested size scaling form (0J. We 
found a good data collapse for all values of p that we 
checked {p = 0.05,.. .,0.5). In Fig. |21 we show the data 
collapse for p = 0.05. This particular value results in the 
best data collapse for the achieveable system sizes. 

We estimated the saturation roughness Wsat ~ 
lim/f^oo w{H) from the flat tail of the roughness curves 
w{H) and show them in Fig.|31as a function of L for sev- 



IV. ROUGHENING IN 3D 

In this section we present our numerical results and 
the corresponding scaling laws obtained for the 3d case. 
If the observations we made in the last section could be 
carried over to the 3d case, one would expect 2 regimes: 
One for small heights H, in which the transverse fluctua- 
tions of the lines are still much smaller than the average 
line-to- line distance a = l/p^^^ in 3d; and one for large 
H, in which the line roughness is restricted due to collec- 
tive effects. Only in the case our line system would also 
for 3d fall into the same universality class as a 3d elastic 
medium (this is, as we have shown the case in 2d), one 
would expect Wsat oc VlnL |^ |^ . 

However, surprisingly we find i) three regimes instead 
of two (c.f. the data shown in Fig. [SJ, and ii) Wgat L, 
i.e. the size of the transverse fluctuations is not restricted 
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FIG. 5: Data for the roughness w; as a function of the height 
H for different transverse system sizes L and two densities: 
p = 0.005 (top) and p = 0.4 (bottom). In the low density limit 
(top) the crossover from single line to collective behavior is 
visible - indicated by the two straight lines with slope 0.625, 
the sinle line roughness exponent, and 0.5, respectively. In the 
high density limit the crossover from collective line behavior 
(indicated by the straight line with slope 0.5) to the saturation 
regime is visible. The data points are averaged over 100-1000 
samples. 




FIG. 6: Scaling plot of the roughness in 3d according to the 
finite size scaling form JSJ for different line densities p = 0.1 
(filled symbols), 0.2 (bold) and 0.4 (empty) for system sizes 
L — 16 (ovals), 24 (squares) and 32 (diamonds). The inset 
shows data collapse according to the scaling form 
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FIG. 7: Saturation roughness as a function of the system 
width with linear fits to the data points. 



by the presence of a large number of other lines but only 
by the lateral system size. Apparently in 3d the lines 
become transparent to each other, and the wandering 
of any line to the transverse direction does not induce 
collective behavior. 

The three regimes that we find can be characterized as 
follows: 1) A single line regime for H <C a^^^ in which the 
roughness behaves as in the the one-line case: w oc 
2) An intermediate regime for a^/^ <C ^ in which 
the roughness increases as H^^^, which is identical to 
the behavior of random walks. Between the two regimes 
one can see a cross-over that can be shown to be related 
to the entropic repulsion of the lines. Recall that this 
leads in 2d asymptotically to collective effects, but here 



the consequences are different. 3) The saturation regime 
for H ^ 'm which the roughness saturates at the 
lateral system size: w k L. In the following we support 
this central result with the data we obtained from our 
ground state calculations for the 3d system and derive 
the appropriate scaling forms for the different regimes. 

In Fig. I^we show our results for the roughness in 3d in 
the crossover region from the intermediate or multi-line 
regime to the saturation regime. We show data for three 
different line density values, but we have also data for 
other values, and they all fit well into the scenario that 
we propose now). The finite size scaling plots yield an 
excellent data collapse using the scaling form: 

w = L-g\^]{H/L^). (6) 
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FIG. 8: Scaling plot for the roughness in 3d in the crossover 
region from single to multi-line behavior. The system size is 
L = 32. The inset shows the original, unsealed data. 



The scaling function g^J{x), which still depends on a, or 
the line density p = has the following asymptotic 

behavior: g\^^{x) cx x^/^ for a; — > and g\^^{x) = const. 
for a; — > oo. 

The first crucial observation here - and the essential 
difference to the 2d case - is that in the limit L — > oo the 
roughness is not significantly restricted by the presence of 
the other lines but approaches a value proportional to the 
lateral system size. Actually, as we see from the plot of 
the saturation roughness as a function of L shown in Fig. 
[7|that Wsat = lim^f^oo w{L, a) — 0.28 ■ L + Ca, where Ca 
is a small constant that varies only slightly with a. This 
variation with a is a boundary effect: The free bound- 
ary conditions act effectively in a repulsive way on the 
lines that competes with the steric inter-line repulsion. 
Therefore systems with a lower line density show smaller 
transverse line fluctuations than those with a higher den- 
sity. 

The second crucial observation is that the roughness of 
the lines in the intermediate regime grows like i.e. 
they have a roughness exponent that is smaller than the 
single line value oi — 0.625 and is identical to the value 
for simple random walks. Although the actual line config- 
uration is constructed in a highly non-trivial manner via 
a global criterion, namely the computation of the global 
iV-line ground state, their universal geometric properties 
appear to be similar to that of random walks. 

The density dependence of the scaling functions g'"^^ {x) 
can be worked out by matching it with the scaling form 
for the single- to multi-line regime. Here the relevant 
length scale in the 7?-direction is a^/'», and in analogy 
to the 2d case we expect for L ^ a the L-independent 
scaling form 

w = a-~gM{Hla^'^) (7) 
with the asymptotics gsdix) = x'' for a; <C 1 and 



g3d{x) = x^l"^ for X \. In Fig. |S1 we show the cor- 
responding scaling plot for the data that we obtained 
from our calculations. Hence we get the expected single 
line behavior w ^ H'^ when H a^/'», and we obtain 

w = a^-(i/2f i7i/2 foj. -u}^L and H > a^/^ . 

(8) 

From this the natural scaling variable for the crossover 
region from the intermediate regime (where w should 
be described by ©) to the saturation regime (where w 
should be proportional to L according to lO) appears to 
be ai-(i/2C)i?i/2/2. or H/a^/^-'^L'^ , which implies that 
© can be rewritten as 

w^L- gsdiH/a^^'^-^L^) for H > a^/^ . (9) 

with gsdix) — x^l"^ for a; — > and gsdi^) — 0.28 for 
X — > oo (see the inset in Fig. |HJ). As in 2d, for high line 
densities p > 0.1 (a < 3) one has to take into account 
the limiting case p — I where lines fill all parallel lat- 
tice bonds resulting zero roughness. This limit can be 
incorporated into lO by rescaling H hy 1/(1 — p). 

At this point we would like to stress that the random 
walk like scaling is not related to the actual distance at 
which the lines touch or cross each other (and are in some 
cases continued randomly) . Both in 2d and in 3d the typ- 
ical length scale s between two consecutive intersection 
points on one line is much larger then the length scale 
f for the crossover from single line to collective behav- 
ior and its divergence with the line average distance a is 
much stronger. This can be seen in Fig. |^ where we 
show scaling plots for the average length of line segments 
s between two crossings, from which one conludes that 
s scales with a as a^'^ in 2d (compared to a-^'^ for the 
crossover lenght scale ^) and a^'^ in 3d (compared to a^-^ 
for ^). This also demonstrates that the lines do not be- 
have like independent random walkers in 3d but reflect 
the effect of a steric repulsion that tends to avoid random 
crossings between them (visible in 2d and in 3d). 

The main result of calculations for the 3d system is 
that the lines with only hard core repulsion can trans- 
verse the whole system, in marked difference to the 2d 
case. One can visualize this result also by looking at the 
disorder averaged position of the center of mass of the 
individual lines (which is r as defined under (PJ). In 2d 
they constitute a regular array on the base line with lat- 
tice spacing a. In 3d, as we show in Fig. ^1 the still 
constitute a regular array, but it concentrates, with in- 
creasing height, more and more in the central region of 
the basal plane of the system. In the limit H ^ oo the 
average center of mass position of each individual line 
will be exactly at the center of the system since nothing 
restricts it from transversing the system from one side to 
the other. 
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FIG. 9: Scaling plots for the length of line segments s between 
two crossings (a) in 2d, the system size is L = 256 and (b) in 
3d, the system size is L = 32. The insets show the original, 
unsealed data. 



WEAK AND ANISOTROPIC DISORDER 
DISTRIBUTION 



Next we check the robustness of the universal scaling 
forms in Eqs. Q, 10, and ^ against a more general dis- 
order distribution. The random energy is distributed 
in [0, S] with 5 > if the bond i is in the transverse direc- 
tion while in [1 — e, 1 -I- e] with < e < 1 if the bond is in 
the longitudinal direction. The disorder distributions on 
the transverse bonds and the longitudinal ones are thus 
different unless e = 1 and S = 2. Moreover, as e becomes 
smaller, the energies on the longitudinal bonds are more 
highly concentrated around 1 implying weaker disorder. 



i'if 



FIG. 10: Average positions of center of masses of lines with 
fixed starting points. From left: H=10, H=100, H=1000. 
Note that the center of mass position is averaged over the 
disorder, here 1000 samples. 



When e = or S^^ — 0, the lines are flat, i.e., w — 
because there is no energy reduction that compensates for 
the energy cost accompanying a transverse fluctuation. 
On the other hand, when e — 1 and S = 2, the lines are 
rough and exhibits the universal scaling behavior shown 
previously. Therefore it would be desirable to find the 
boundary between such distinct two limits in the (e, S)- 
plane. We computed the roughness for the generalized 
disorder distribution with various values of e and (5, which 
is shown in Figs. (2d) and 1121 (3d). We consider for 
2d systems with the fixed line density N/L = 1/16. In 
3d, a single line is studied to sort out the effect of e and 
6 on the roughness scaling in Eqs. Q), iQ, and 10 more 
clearly. The roughness shows the same scaling behaviors 
as in the previous sections except for the dependence on e 
and S before it saturates. It implies that interacting lines 
are rough for nonzero e and and that the scaling 
behaviors in Eqs. (QJ, ((SJl, and © are universal with 
respect to the variations of the disorder distribution used 
here. 

The dependence of the roughness on e and 6 before 
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FIG. 11: Data collapse of the scaled roughness w/lnL vs. 
the scaling variable H(e/S)'^ /(In L)^^'' of interacting lines in 
the 2d lattice of size L x H with L = 16, 32, 64, and 128 
and the line density N/L — 1/16. The parameters e and S 
are chosen as in the legend. The exact wandering exponent 
= 2/3 |]3 and the fitting value v ~ 0.8(1) are used for the 
scaling variable. The inset is the plot of the roughness w vs. 
H for L — 64. Dependence of w; on e and 5 is shown before 
the roughness saturates. 
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crossover scaling function between the transient regime 
and the stationary state is represented in terms of the 
scahng variable H / / (wsa.t{L)Y/'^ , which is shown in 
Figs. ^2 and E| with the fitting value for v in ~ 
{t/6y being 0.8(1) (two dimension) and 0.68(8) (three 
dimension), the latter being in rather good agreement 
with 2/3. The derivation of (|10|1 corresponds to compar- 
ing elastic energy to the pinning energy in a conventional 
derivation of the Larkin length. Obviously in two dimen- 
sions one should be more careful while simply equating 
the costs of tranverse bonds to the elastic energy of the 
line and the costs of longitudinal bonds to the strength 
of a pinning potential. 



H (e/5)7L 



FIG. 12: Data collapse of the scaled roughness w/L vs. the 
scaling variable H{e/&Y / L?'^'' of one line in the 3d lattice of 
size L X L X H with L = 8, 16, and 32. The wandering 
exponent C = 0.625 [H d is used and the exponent u is 
found to be 0.68(8). The inset is the plot of the roughness w 
vs. H ior 8. 

it saturates is rooted in the crossover behavior around 
the point (e,(5~"'^) = (0,0). As {e,S~^) deviates from 
(0,0), a line roughens {w > 0), if the reduction of the 
longitudinal-bond energy obtained by the transverse fluc- 
tuation is larger than the accompanying energy cost. By 
comparing the energy reduction and cost accompany- 
ing the transverse fluctuation, one can get the crossover 
("Larkin") length scale for H, For a given line i, 
the typical fluctuation or the standard deviation of the 
sum of the energies on the longitudinal bonds, E'^ || = 

J2z e||(ri(2:)), is given by e^/ll, which corresponds to the 
average reduction of || achieved by a transverse fluc- 
tuation. On the other hand, the line chooses the posi- 
tion where the cost is minimum among H possible ones 
(1 < z < H) to fluctuate in the transverse direction. 
The cost on the average is equal to the average mini- 
mum value of H independent random variables each of 
which is distributed uniformly in [0, S], which is given by 
S/{H+1) [1^. Comparing the energy reduction and cost, 
one finds the characteristic length scale as 



such that the roughness w is zero for H ^ and non- 
zero, increasing as H'^ for H ^ £^h- From the numerical 
data and the above argument, the roughness in 2d and 
3d for the generalized disorder distribution used here can 
be represented as 

{§y ii«§«ML)]'/<), (11) 

[ WsAL) (^»Kat(i)]^/0, 

with the usual roughness exponents and saturation 
roughness scalings according to the dimension. The 



VI. CONCLUSIONS 

In this work we have considered the collective effects 
arising from the interaction of "forests" of directed poly- 
mers, in the presence of quenched randomness. In both 
2d and 3d this is marked by a cross-over from individ- 
ual, rough lines to a "collective" regime with a system- 
size dependent roughness w. In 2d, our studies augment 
earlier work on similar systems using fixed line densi- 
ties. The approach is based on a mapping of the prob- 
lem to the minimal matching problem of combinatorial 
optimization, and the essential point of our model(s) is 
that the density of lines is a free parameter. In this, 
we consider the cross-over to the asymptotic collective 
roughness as a function of the parallel and perpendicular 
disorder strengths. 

In 3d, perhaps more importantly, it is found that 
the generally expected logarithmic roughening is absent: 
while (scalar) models of weakly perturbed elastic mani- 
folds do exhibit logarithmic roughness scaling this is not 
the case here. An ensemble of directed polymers with 
hard-core interactions exhibits almost "trivial" roughen- 
ing, with a random-walk -like scaling picture. There is 
a crucial difference between two and three spatial di- 
mensions even within the confines of the same micro- 
scopic model. An ensemble of directed polymers with lo- 
cally repulsive interactions is thus proven to be different 
from other approaches using mappings to e.g. higher- 
dimensional manifolds that do exhibit logarithmic corre- 
lations 9]. It is worth noting that this does not follow 
from a "profileration of dislocations" in a vortex lattice, 
but from the absence of collective behavior. 

A rather random but likewise potentially very impor- 
tant consequence of the same is that in 3d the lines en- 
tangle - the topological state becomes highly non-trivial 
such that it might be useful to describe it in terms of 
knot-theoretical means 0, 0. Such concomitant ge- 
ometrical structures are again absent in models for 3d 
elastic media. In this case, the description of barriers 
and excitations still remains to be done including the ef- 
fect of the entanglement on both. 
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